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Abstract 

Nonclassicality in the sense of quantum optics is a prerequisite for entanglement in multi-mode 
radiation states. In this work we bring out the possibilities of passing from the former to the 
latter, via action of classicality preserving systems like beamsplitters, in a transparent manner. 
For single mode states, a complete description of nonclassicality is available via the classical theory 
of moments, as a set of necessary and sufficient conditions on the photon number distribution. We 
show that when the mode is coupled to an ancilla in any coherent state, and the system is then 
acted upon by a beamsplitter, these conditions turn exactly into signatures of NPT entanglement 
of the output state. Since the classical moment problem does not generalize to two or more modes, 
we turn in these cases to other familiar sufficient but not necessary conditions for nonclassicality, 
namely the Mandel parameter criterion and its extensions. We generalize the Mandel matrix 
from one-mode states to the two-mode situation, leading to a natural classification of states with 
varying levels of nonclassicality. For two-mode states we present a single test that can, if successful, 
simultaneously show nonclassicality as well as NPT entanglement. We also develop a test for NPT 
entanglement after beamsplitter action on a nonclassical state, tracing carefully the way in which 
it goes beyond the Mandel nonclassicality test. The result of three-mode beamsplitter action after 
coupling to an ancilla in the ground state is treated in the same spirit. The concept of genuine 
tripartite entanglement, and scalar measures of nonclassicality at the Mandel level for two- mode 
systems, are discussed. Numerous examples illustrating all these concepts are presented. 

PACS numbers: 03.67.Bg, 42.50.Ar, 42.50.Ex, 03.67.Mn, 03.67.-a 
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I. INTRODUCTION 



States of multi-mode quantized radiation fields are, by very definition, nonclassical in 
nature. However for many purposes it is useful to characterize certain states as displaying 
nonclassical features in a particularly prominent or manifest manner. The quantum optical 
concept of nonclassicality, based on the diagonal coherent state representation [1], and the 
notion of entanglement [2], are two such important nonclassical features displayed by some 
states but not by others. Whereas the former has played an important role from the early 
days of quantum optics [1, 3-12], the latter has received enormous attention more recently, 
with the development of the theory of quantum information [13-17]. Of course entanglement 
is meaningful only in the case of two or more modes of radiation, while nonclassicality is a 
useful concept even at the single-mode level. 

The main aim of this work is to study the relationships between quantum optical non- 
classicality and entanglement of various kinds, and the possibility of generating the latter 
given the former. In this context the asymmetry between these two features must be care- 
fully appreciated. It is well known that quantum optical nonclassicality is a prerequisite for 
entanglement [18]. Whereas every entangled state is nonclassical, a nonclassical state may 
be separable (tensor product of nonclassical states, for instance) or entangled. Thus entan- 
gled states are a proper subset of nonclassical ones. Then the following questions become 
meaningful. For every signature of nonclassicality, how much further must one go and what 
additional conditions have to be met before one is assured of entanglement (of specified 
type)? Given a nonclassical separable state, can it be transformed into an entangled state 
through a physically realizable process using classicality preserving systems? 

In speaking in this sense of converting nonclassicality into entanglement we have in mind 
the use of passive devices such as beamsplitters [19, 20, 34] which act on several modes 
of a given system to produce unitary linear combinations of the annihilation operators of 
the modes. Such systems preserve total photon number and, further, they cannot create 
or destroy nonclassicality. A general study must include developing and clearly specifying 
signatures of various levels of nonclassicality on the one hand, and of entanglement on the 
other, and then examining how one can reach the latter starting from the former, using 
these classicality preserving passive systems. 

Two comments by way of clarification are in order before we describe the manner in 
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which the contents of this paper are organized. First, it should be stressed that our concern 
in this paper is (almost exclusively) with non-Gaussian states. Nonclassicality of Gaussian 
states is necessarily of the squeezing type (quadrature squeezing), and a beam splitter can 
convert such a nonclassicality into entanglement. It is essentially this Gaussian scenario 
that largely formed the work-bench during the early years of the theory of entanglement 
for continuous variable (canonical) systems, elevating the role of the symplectic group of 
linear canonical transformations and its unitary (metaplectic) representation. The interplay 
between nonclassicality and entanglement in the Gaussian case has been long appreciated. 
For instance, the principal result of Ref. [15] can be phrased thus, as emphasised by Kimble 
recently [21]: A two-mode (mixed) Gaussian state is separable if and only if its nonclassicality 
can be removed by local linear canonical transformations — entanglement is the same as 
nonlocal nonclassicality. 

More recently, however, there has emerged considerable interest in non-Gaussian states 
[22-33]. Ref. [31-33] may be consulted for a good introduction to the literature. Nonclas- 
sicality signature of non-Gaussian states is not restricted to quadrature squeezing — the 
state can be antibunched, for instance. A non-Gaussian state which is neither squeezed nor 
antibunched could exhibit nonclassicality signature in the higher moments of the photon 
count statistics or photon number distribution [9]. Non-Gaussian states thus exhibit a rich 
variety of nonclassicality signatures. 

The second comment is in respect of beam splitters; we make extensive use of this gadget 
in the present work. In the domain of quantum optics, nonclassicality generating operations 
like squeezing are considered to be expensive resources whereas beam splitter, being a passive 
device, is rightly considered inexpensive. The latter is passive in the sense that it cannot alter 
the total number of photons in the pair of modes it couples, and also in the stronger sense 
that it cannot create or destroy nonclassicality. It thus makes sense to ask to what extent 
the expensive nonclassicality resource can be converted into the (expensive) entanglement 
resource using the inexpensive beam splitter resource. One should, however, appreciate that 
the beamsplitter effects a joint unitary transformation on the pair of modes it couples, and 
hence will be considered an expensive resource in a different context like a pair of nano 
mechanical oscillators. Our analysis thus applies to the specific context of quantum optics. 

The contents of this paper are arranged as follows. Section II defines the concepts of 
quantum optical classicality — QO-cl and nonclassicality — QO-noncl, for general single- 
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mode fields. In the phase invariant case they are entirely stated in terms of the photon 
number distribution [35]. This is possible thanks to the result of the classical Stieltjes 
moment problem [36]. It is then shown that if such a single- mode state is coupled to an 
ancilla in any coherent state and passed through any nontrivial U(2) beamsplitter, the 
resulting two-mode output state shows NPT entanglement [37] precisely when the input 
single-mode state is quantum optically nonclassical [38, 39]. In this case the signatures of 
the two coincide exactly. Section III defines quantum optical classicality and nonclassicality 
for general two-mode states. We then describe a single test which if successful is able to 
establish simultaneously both the nonclassicality and NPT entanglement of a given two- 
mode state. Of course proof of entanglement automatically implies proof of nonclassicality, 
so the interest here is in the structure of the expressions used in the test. Since theorems 
of the moment type are not available for two or more modes, we describe in Section IV a 
sufficient criterion for one- mode nonclassicality originally due to Mandel [8] , and extend it to 
two-mode systems as well. All the later considerations of this work are based on this Mandel 
level of nonclassicality which is of course weaker than the general notion of nonclassicality. 
Several levels or types of nonclassicality in this sense are described, and three examples to 
illustrate the ideas are presented. 

In Section V we study states of two-mode systems for which both notions of nonclassicality 
and entanglement make sense. We find conditions for nonclassical states to become NPT 
entangled via BS action, and show in detail the way in which these conditions go beyond 
those that guarantee Mandel type nonclassicality. This analysis is applied to two examples 
to see the formalism in action. Distillability of the resulting state is demonstrated in one 
case. Section VI takes up the theme of Section II and extends it to higher number of modes. 
Thus we couple a two-mode state which is nonclassical to a third ancilla mode in vacuum, 
pass this three-mode state through a beamsplitter acting on the three modes, and develop 
tests for NPT entanglement in the output state. As in Section V, here too the precise manner 
in which testing for entanglement goes beyond testing for nonclassicality of the Mandel type 
is emphasized. Section VII discusses the subtle notion of genuine tripartite entanglement, 
while Section VIII brings out some features of two-mode Mandel level nonclassicality and 
beamsplitter action. All the ideas in Sections VI, VII, and VIII are illustrated through 
several examples with tractable analytical structures. Section IX contains some concluding 
remarks. 
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II. SINGLE-MODE FIELDS 



Let the photon annihilation and creation operators for the concerned mode be written as 
a and a* , obeying the commutation relation 

[a,a f ] = 1. (2.1) 

The familiar Fock states and coherent states are 

\n) = (n!r 1/2 (at)», „ = 0,1,2, ••• ; 

1 2 °° Z Ha 

\z a ) = e-*l*«l Y,-f=j\n a ), a\z a ) = z a \z a ), z a eC. (2.2) 

n a =0 V n a- 

The Fock states {\n)} form an orthonormal basis for the space of all single-mode states, 
while the coherent states {|z a )} form a normalized nonorthogonal overcomplete system. 

For a single mode there is no meaning to entanglement, only a separation of states into the 
quantum optical classical (QO-cl) and the quantum optical nonclassical (QO-noncl) types. 
This is based on the diagonal coherent state representation of a general (pure or mixed) 
state p (a) : 

pW^Tr- 1 / d 2 z a <f>(z a )\z a )(z a \. (2.3) 
J c 

The properties of the real diagonal representation weight 4>(z a ) determine the nature of p^ 
[1]: 

4>{z a ) > & p (a) QO-cl; 

4>{z a ) 2 p {a) QO-noncl. (2.4) 

In the former case, the meaning is that (f>(z a ) is pointwise nonnegative; it is then mathemat- 
ically a valid probability density in phase space. 

If one is interested only in the expectation values of number conserving observables, i.e., 
of operators commuting with N a = a^a, and hence diagonal in the Fock basis, it suffices to 
use the phase averaged form p$ of p^ : 



4a) 

Pd 



e iBN a p{a) e -i9N a 



Jo 2tt ( 
= - d 2 z a P{I a )\z a )(z a \, I a 

71 J 
r 2w Aft 



NaP^- 



(2.5) 
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Clearly is a physical state, equivalent to p^ as far as expectation values of number con- 
serving operators are concerned. Moreover, the quantity P(I a ) is the diagonal representation 
weight for pffl, so (2.4) leads to the coarse grained classification [11, 12] 

P(Q > & p [ n ] QO-cl; 

P(Q t p { $ QO-noncl. (2.6) 

Further, all information about P{I a ) is contained in the photon number probabilities p(n) = 
(n\p^\n) = (n\p$\n) for n — 0, 1, 2, • • • [35]. It is convenient to call the set of probabilities 
{p(n)} the photon number distribution (PND). Clearly, p^ 1 can be explicitly written in 
terms of the PND and vice versa : 

oo 

p ( d = Y,p( n )\ n )( n \> 

n=0 

p(n) = / dI a P(Qe-^I a n /n\ > 0. (2.7) 
Jo 

Therefore one can ask if the conditions (2.6) for the coarse grained QO-cl — QO-noncl 
divide can be explicitly given in terms of the PND {p(n)}. This is indeed possible [35], as 
a result of the classical analysis of the Stieltjes moment problem. It involves two infinite 
sequences of matrix positivity conditions, set up as follows. It is convenient to introduce the 
auxiliary quantities 

q n = n\p(n), n = 0,1,2,- (2.8) 
Then define two real symmetric infinite dimensional matrices L and L as follows : 

L = (L n ' n ), L n i n = q n '+ni n ,n = 0,1,2, ■■■ ; 

L = (L n / n ), L n > n — q n / +n+1 , n',n = 0,1,2,- ■■ . (2.9) 

The diagonal elements L nn of L are q 2n , while L nn of L are (fcn+i- Using L and L, we can 
define a sequence of (N + l)-dimensional matrices L^ N \ as restrictions. Then the key 
result is [35] : 

P(Q > {p { d ] and PND {p(n)} are QO-cl} 
& L {N \ P>0, JV = 0,1,2,...; 
i (JV) = {Ln'n : n',n = 0,l,2,---,N), 

LW = (L n , n : n',n = 0,l,2,---,N). (2.10) 
7 



Both L( N 1 and are real symmetric (N + l)-dimensional matrices, made up of the 

intersections of the first (N+l) rows and (N+l) columns of L and L respectively. Conversely 
we have : 

P(I a ) £ <£> and PND {p(n)} are QO-noncl} 

<^> some L {N) £ and or some L {N) £ 0. (2.11) 

As noted earlier, for one-mode states there is no concept of entanglement. However a QO- 
noncl state characterized by (2.11) can be converted to a state possessing NPT entanglement 
by adjoining it to a second ancilla mode, the 6-mode, initially say in its vacuum state, and 
passing this two-mode state through a generic beamsplitter (BS) corresponding to an element 
of the group U{2). This is seen as follows : 

The 6-mode brings in an added operator pair b, obeying the same commutation relation 
(2.1); and each of a, a) commutes with b and b'. Let now 

u = 

\U 2 1 U 22 

be a general element of the two-dimensional unitary group U(2). The action of the cor- 
responding BS on the operators of the two modes is given by conjugation with a unitary 
operator U acting on the two-mode Hilbert space, the dependence of U on u being left 
implicit [11, 34] : 

U" 1 I " I V = u f " I , U~ l f " I U = u* [ a I . (2.13) 
{„) \i) {») {»} 

This unitary operator is 'passive' in the sense that 

U(N a + N b ) = (N a + N b )U. (2.14) 

It is also passive in another important sense in the present context : for any two-mode state 
the property of being QO-cl or QO-noncl is preserved upon passage through any BS, so BS 
action can neither create nor destroy nonclassicality. ( The output state of a beam splitter 
is a coherent state if and only if the input is. This will be clarified and elaborated in the 
next Section). 
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(2.12) 



We start with the a-mode state p^ which commutes with N a , and take as the two- mode 
input to the BS the separable (product) state 

p! n a6) =#®|0) b& (0|. (2.15) 

This commutes with the total number operator N a + iVj,. Passage through the BS preserves 
this property and results in the output state 



= llf:^^T\0,0)(0,0\(arU- 



n=0 nl 

00 p(n^ 



= E^r( M n« f + ^1^)10, 0)(0,0|K 1 a + ^ 1 6) n . (2.16) 



n=0 nl 



As a general notation let us now use n, n', n", ■ ■ ■ to denote the eigenvalues of N a , and 
to, m',m", ■ ■ ■ to denote eigenvalues of 7V b . Then the general Fock state matrix elements of 

4ab) 

Pout are • 

n' m' *n *m 

I 1 i\4ab)\ \ r u ll u, 21 "ll u 21 /r> ,71 

(n ,to |p out J 77,, to) = ^ +m>+m g n+m —==—==. (2.17) 

Vfi !m ! vto!to 

To see whether p^ possesses NPT entanglement, we carry out the partial transpose (PT) 
operation by implementing it in the 6-mode space in the Fock basis. This amounts to 
interchanging to and to' on the right hand side of Eq. (2.17), and gives: 

„,n'„,ra „,*n„,*m' 

An' rr,'\n^ PT \r, rn\ — A n M 11 M 21 M 11 M 21 Icy ■> n\ 

[n , m |p out | jt,, to) - <V+m,n+m' g n+m / —===—====. (^• i °) 

\m ™! vto'to' 

The question now is whether this is a physical state i.e., if p Q u? PT > 

We now isolate two principal submatrices out of the matrix p^ PT in the Fock basis, 
which are closely related to the matrices L, L of Eq. (2.9). The first one is the submatrix 
H = (H n r n ) obtained by taking m' = n', m = n in (2.18) : 

H n , n = {n,n \p y out J | n, n) 

(u* u u 21 ) n 

= Qn+n* y. j , n,n = 0, 1,2, 

n'l nl 

H = A^LA, A = diag( ( ^ lM , 2l)n , n = 0, 1, 2, • • •). (2.19) 

71! 

The second is the submatrix H = (H n ' n ) obtained by taking to' = n' + 1, to = n + 1 in 
(2.18): 

£r ///ill 4ab)PT } . \ 

ii n / n = (n,n +l|p ou t |n,ra + l) 



{u u u* 21 ) n ' u* 21 (u* llU21 ) n 



u 



gn+n'+i 7i , = , = ; ^=F, n,n = 0,1,2, 



n'! 



vV + 1 n! V« + 1 ' 



H = A^LA, A = diag( ( "* 1 " 2l)W n = 0, 1, 2, • • •)• (2.20) 

n! Vn + 1 

Invertibility of A and A implies that H, H > if and only if L, L > 0. We then see: 
every signature of and the PND {p(n)} being QO-noncl, such as some ^ or 

some L( N ~) 0, directly implies a corresponding signature of p^ut being NPT entangled, 
since either H ^ or H ^ 0, implying that p^t PT ^ 0. It is in this precise sense that 
any nontrivial BS is able to convert input QO-nonclassicality of a single-mode to NPT 
entanglement of the output state p^ , in a sense preserving the signature of nonclassicality. 

A simple calculation shows that in the input state (2.15) of the a6-system, the 6-mode 
need not be in the vacuum but could be in a general coherent state \z b ) b for some complex 
nonzero z b . Since 

\z b ) b = D^(z b )\0) b , 
D^(z b ) = exp(z b tf -zfb), (2.21) 

we have the following replacements for the previous Eqs. (2.15, 2.16) : 

Pin &) = P ( D ® \z b )bb{z b \ 

= D^ b \z b ){p { S ® |0> w (0|}DW(^) _1 ; 

pS = Up^u- 1 

= D (a \u 12 z b )D^ b \u 22 z b ){Up^ ® ^^(OlU-^D^^z^D^iu^y 1 . (2.22) 

(In passing we note that this input state no longer commutes with N a + N b ). Here we 
have used Eq. (2.13). Thus the matrix elements of this two-mode BS output state in the 
displaced orthonormal product basis D^ a \ui2Z b )D^(u22Z b )\n ) m) are identically the same as 
the matrix elements of of Eq. (2.16) in the Fock basis \n,m). Implementing the PT 
operation in this new product basis we recover the earlier results. Thus every signature of 
QO-noncl of p$ is transformed by a general BS action into a corresponding signature of 
NPT entanglement of the output two-mode state (2.22). 

In both cases (2.16) and (2.22) we see that nontrivial BS action, while being passive 
and so maintaining any QO nonclassicality present in the input state, is able to convert an 
initially unentangled two-mode state into an NPT entangled state. This it can do provided 
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there is some QO nonclassicality to begin with. Crudely speaking, what the beam splitter 
achieves in the present case is to take the input nonclassicality which resides 'locally' in the 
a-mode and convert it into nonclassicality residing 'nonlocally' as entanglement between the 
modes. 

III. TWO-MODE FIELDS - GENERAL PROPERTIES AND AN ENTANGLE- 
MENT TEST 

We now consider specific new features encountered in the study of states of a two-mode 
system. A general two-mode state p^ possesses the diagonal coherent state representation 

= I lt^t^^ Zh) \ z ^ Zh){z ^ Zh \ (3.!) 

J J 71 71 

in terms of the two-mode (product) coherent states \z a ,z b ). Analogous to Eq. (2.4), the 
properties of (f>(z a , Zb) determine the nature of p^ : 

<f>{z a , Zb) > & p {ab) QO-cl; 

(f)(z a , z b ) t <=> p {ab) QO-noncl. (3.2) 

In the former case <f>(z a , Zb) has the properties of a probability distribution in the four 
dimensional phase space. 

Now, apart from examining whether a given state p^ is QO-cl or QO-noncl, we can also 
ask in the latter case whether it is entangled, and if so whether it is NPT type, distillable [40] , 
etc. These added questions become meaningful. In fact we will develop later in this Section 
an interesting test or criterion which can witness simultaneously for QO-noncl assicality as 
well as NPT entanglement. 

With respect to BS action (2.12, 2.13) representing general elements u G U(2), we note 
the following. Such action is nonlocal since the modes a and b get linearly mixed, in addi- 
tion to being passive in the sense of conserving N a + N b . Since annihilation operators go 
into linear combinations of annihilation operators under this action, coherent states go into 
coherent states with the matrix u acting on the pair z a , Zb as a column vector. Therefore 
convex sums of coherent states go into similar convex sums [41], (j>(z a , Zb) experiences a point 
transformation <f>(z a , Zb) — > <f>(uuz a + ui2Zb, u 2 iz a + U22Zb), and thus such BS action preserves 
the QO-cl or QO-noncl nature of the state p( ab \ This is another sense of BS action being 
passive. 
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On the other hand, while a QO-cl state has no entanglement, a QO-noncl state may be 
separable, i.e., unentangled, or may possess entanglement. Entangled states are a proper 
subset of QO-noncl states, and NPT entangled states are a further subset. Now BS action 
can cause a transition, within the QO-noncl subset, from a separable to an entangled state, in 
which case we can further enquire into the nature of the entanglement so obtained, whether 
it is NPT type etc. This is in fact the case in the transition from p[^ of Eq. (2.15) to p^t 
of Eq. (2.16), or the transition in Eq. (2.22). 

Continuing with the definitions of QO-cl and QO-noncl in Eq. (3.2), if we are interested 
only in the total number conserving matrix elements of various operators, i.e., of operators 
commuting with N a + N b , it suffices to work with the total phase averaged state 

Jab) f 27T ^_i0(N a +N b ) Jab) -iO(N a +N b ) 

Pd ~ Jo 2n e 9 6 

d 2 z a d 2 z b 



= f _?l_?±p(j aj I b ,8)\z a ,z b )(z a ,z b \, 

J 71 71 



h = \z a \ , h = \z b \ , 9 = argz* a z b , 
P(h, h, 9) = r ^(v£e-*\ v^-*')). (3.3) 

This state is number conserving : 

p ( f(N a + N b ) = (N a + N b )p% b \ 

(n'm'\p^ b) \nm) = 5 n > +m ^ n+m (nm\p {ab) \nm). (3.4) 

Since the coarse grained P(I a , h, 9) is the (real) diagonal representation weight of pfj , we 
have the following QO classification at this level [11, 12] : 

P(I a ,h,0)>0 & P ( D b) is QO-cl, 

P(I a , h, 9) ~£ & p ( £ b) is QO-noncl. (3.5) 

We use this in the next Section. 

We now revert to a general state p^ and describe a test which, if it succeeds, simulta- 
neously establishes both QO nonclassicality of p^ and its NPT entanglement. It is clear 
of course, that any test which establishes the latter immediately implies the former; the 
interest here is in the structure of the test itself. 

We set up an infinite matrix N with operator entries Njk,i m where j, k, I, m run indepen- 
dently over 0, 1, 2, • • •. The pair jk denotes a 'row index' and takes in sequence the values 
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00; 10, 01; 20, 11, 02; 30, 21, 12, 03; • • • . Similarly the 'column index' pair Im also takes these 
same values in the same sequence. We define the entries of N thus : 

N jkM = JV,JU = a«St*a'S™. (3.6) 

Clearly N = (Njk,irn) is an infinite 'hermitian' matrix of operator entries. Note that 
these entries are in normal-ordered form. Starting with the diagonal representation (3.1), 
for any set of complex coefficients {cj k } and the associated positive semidefinite operator 
Y,jk,im c jkNjkM c irn, we always have: 

Tr(p^ J2 c* jk N jkM c lm ) = Tr(p^(^c, fc a^ fe )t(^ Qm a / 6 m )) 

jk,lm jk Ira 

= I / — — IE > 0. (3.7) 

J J TV 7T , 

Ira 

This is independent of (f>(z a , z^) being classical or otherwise, because we have here the ex- 
pectation value of a positive semidefinite hermitian operator. On the other hand, if we pass 
to the partial transpose p( a6 ) PT of p^, by performing transposition only in the space of 
states of the 6-mode in the Fock basis, this will amount to everywhere replacing b^b m by 
b^ti , since in the Fock basis & and b are real [42]. Thus for the same positive semidefinite 
operator as in (3.7) we have: 

TY{p (a b )PT £ c * jk N jkM c lm ) = Tr(p^) £ c* jk a^P m a l b k c lm ) 

jk,lm jk,lm 

dlZb <t>{z a ,z h )\Y,c lm z l a 4 m \ 2 - (3-8) 



// 



7T 7T 



Ira 

Notice the difference in the integrands of the last integrals in (3.7) and (3.8); the latter 
integral is sure to be positive if p( afe ) PT > 0, otherwise it could be negative. 

Thus we arrive at a single step test for QO-nonclassicality and NPT entanglement of 
p(ab) _ a k ove expression (3.8) being negative implies two things simultaneously: 

(i) <f>(z a , z h ) ^ 0, hence p^ is QO — noncl; 

(ii) p {ab)PT £ 0, and hence p {ab) is NPT entangled. (3.9) 

As we said already, p^ being NPT entangled already implies also its being QO-noncl. The 
point here is that the expression (3.8) being negative manifestly displays both properties of 
p( ab ) immediately. 
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This interesting result is an indication of the possibility, in suitable circumstances, of 
'bridging the gap' between the characterization of QO nonclassicality and the further char- 
acterization of (NPT) entanglement for two-mode fields in an efficient manner. For the 
particular the test based on (3.8), there is no 'gap' at all, but as we will show this will not 
always be the case. 



IV. MANDEL MATRICES AND STATE CLASSIFICATION FOR ONE AND 
TWO-MODE FIELDS 

The discussion of single- mode states in Section II made use of the probabilities p{n), 
n = 0, 1, 2, • • • , individually. In terms of the moments of the PND, this means in principle 
that all its moments are involved. It is naturally possible to make much more limited 
statements about QO nonclassicality if one uses only the first and second moments, say, of 
the PND. This is the content of the Mandel criterion for QO nonclassicality stated in terms 
of the Mandel matrix (or Mandel Q parameter) associated with a given state. It is useful 
to begin by outlining this for single-mode states to set up notations, and then generalize to 
two modes. 

Another motivation is that for two or more modes there are no moment theorems at 
all comparable in scope to the one stemming from the classical Stieltjes moment problem, 
and therefore no simple generalizations of Eqs. (2.10, 2.11) into necessary and sufficient 
conditions for QO classicality can be expected to be available. 



A. Mandel matrices for single-mode states 

Given a one-mode state p^ a \ we construct by Eq. (2.5) the state p^ 1 conserving N a . The 
fact that pjj is a physical state leads to a certain 2x2 real symmetric matrix, involving the 
first and second moments of the PND {p(n)}, being positive semidefinite : 

#>>«>=► | 1 = ( 1 <a ' a> |>o, 

(Ay W) ) \ (8*8) (a'asta) 

oo 

(N a ) = Enp(n), <^> = £n 2 p(n). (4.1) 

n=0 n=0 
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All expectation values here are in the state p$ . The Mandel matrix associated with p^ 
is obtained by replacing the expectation value (N 2 ) by the expectation value of the normal 
ordered form a) 2 a 2 of N 2 : 



( 



(a)^ 



(4.2) 



1 (ata) 
y (at a) (at 2 a 2 ) f 

Here the superscript ^ indicates that we are considering states of a single-mode. While 
positivity of p$ implies the positivity of the 2x2 matrix in (4.1), it does not imply the 
positivity of M^(p^). From Eq. (2.7) we easily obtain: 

(at a) = 5>p(n) = / dI a P(I a )I a = (I a ), 

°° f-OG 

(at oat a) = Y, n 2 p(n) = / dI a P(I a )I a (I a + 1) = (I a (I a + 1)), 

°° roo 

(a^a 2 ) = Y,n(n- l)p(n) = jf dI a P(I a )I 2 = (I 2 ), 



n=0 
roo 



(a*a 2 ) - (tfa) 2 = / dI a P(I a )(I a - (I a )) 2 . 
Jo 



(4.3) 



Clearly, the last expression cannot be negative if P(I a ) is pointwise nonnegative. Thus we 
are led to the Mandel classification of states of a single-mode field [8, 12] : 

P(Q > 0, p { $ is QO-cl {M (1) (p^) > <£> detM^(p^) > 

& (AN a ) 2 - (N a ) > 0} : 

Pc? displays super- Poissonian statistics (super-PS); 



{M\p { $) t & detM^(p^) < & (AN a ) 2 - (N a ) < 0} 
P(Q t 0, P ( S is QO-noncl : 

p$ displays sub- Poissonian statistics (sub-PS). (4.4) 

It is because the Mandel matrix is two dimensional with obviously positive trace that the 
positivity or nonpositivity of MW(^') reduces to that of its determinant, hence to that of 
(AN a ) 2 -(N a ). 
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B. Single-mode squeezed vacuum example 



As a useful and instructive example of this concept, we consider the case of a single- mode 
squeezed vacuum [3-5]. Such a state is obtained by applying a unitary (scaling) operator 
involving the exponential of a complex combination of a) 2 and a 2 to the Fock vacuum |0) o , 
and is parametrised by a complex variable £ = (1 + ^2 oran equivalent complex variable u : 

\^ a \u)) = exp{i(£at 2 -ra 2 )}|0) a 



(i - M 2 ) J T. 



n=0 

u> = ||tanh(|£|/2). (4.5) 

Since only even photon number states are present, the probabilities p(l), p(3), p(5), • • • in 
the PND vanish. That is, L ^ but all its diagonals vanish, implying L ^ 0, which is 
immediate evidence that these states are QO-noncl. Some important expectation values 
are : 

(^(u)\{a\a, N a , Nl <Pa 2 , a 2 }\^(co)} = 

{0, 0, S 2 , S 2 (S 2 + 2C 2 ), S 2 (2S 2 + C 2 ), 
S = sinh(|£|/2), C = cosh(|£|/2). (4.6) 

The 2x2 Mandel matrix for this state is thus : 



M«(|V (a) H)) 



1 s 2 \ 



S 2 S 2 (2S 2 + C 2 ) ) 
detM«(|^ (a) (w))) = S 2 (S 2 + C 2 ) > 0, (4.7) 

where S and C are given in Eq. (4.6). Thus these states have super-PS, and the QO- 
nonclassicality does not show up, or is missed, at the Mandel level. 



C. Mandel matrices for two-mode states 



The two-mode generalization of the Mandel matrix idea leads naturally to a more intricate 
classification of states. We consider only states pffi conserving, i.e., commuting with, N a + 
N b . First we develop the analogue of the positivity property (4.1), the statement that the 
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number fluctuation (AN a ) 2 for a PND {p(n)} can never be negative. Define a column and 
row vector with number conserving operator entries as follows 



C 




C 4 = (iV a S+a a+S iV 6 



(4i 



With their help next define a 5 x 5 matrix with operator entries and which is 'hermitian' 
like N in Eq. (3.6), and also 'positive definite' : 

1 ct 



( i 



1 C 1 " 
(7 CC* 



\ 



(4.9) 



TV 2 iV a 6ta iV a at6 A^iVf, 

atft a+SiVo a^a (a f 6) 2 a+SiVi, 

fota 6t ai v a (ftta) 2 fetaatft paN b 

\N b N b N a N b tia N b tfb N b ) 

Given a state fffi, by taking entrywise expectation values of the operators in E we get 
the 5x5 numerical hermitian matrix E which is clearly hermitian positive semidefinite and 
generalizes (4.1) : 

' 1 ct 



E = (E) = Trffi 



(ab) 



( 



(Ct) 



> 0. 



(4.10) 



V 



(C) (CCt) 

We now define the two-mode Mandel matrix for the state by replacing C& in Eq. 
(4.10) by its normal ordered expression (the entries in C and & are already in the normal 
ordered form) [12] : 



B 



{ab). 



: C& :, 
Tr(pif 



( x fr\ 



C B 
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( 1 (ata) (ft a) (atfe) (&t&) ^ 

(ata) (at 2 a 2 ) (^a 2 ) (tf 2 ab) (tfftab) 

(ati) (at 2 a6) (a^afe) (at 2 6 2 ) (atfe+6 2 > 

(ftta) (a^a 2 ) (6t 2 a 2 ) (a^) (&t 2 a&) 

\ <S+S> (atfttaS) (ti 2 db) (tfftP) (bH 2 ) j 

= di a di b —p(i a ,i b ,e) x (i a ,h,o)x(ia,h,oy 

JO JO Jo ATI 

X(la,h,0) = (l I a JTj h e+ i(> VTj b e-^ I b ) T (4.11) 

The superscript < - 2 - ) indicates that we are dealing with a two-mode state, and this Mandel 
matrix is 5 x 5 hermitian but not necessarily positive semidefinite. 

Before presenting a classification of two-mode states based on the 5x5 Mandel matrix, we 
introduce a useful derived object. This is the 2x2 Mandel matrix associated with a general 
single mode defined as a linear combination of the modes a and b, the reduced subsystem 
state of this chosen single mode being calculated with respect to the two- mode state (f^. 
The definition of the annihilation operator A of such a mode and then of its Mandel matrix 
are : 



A = ad + (3b, 
M™(p^;a,(3) 



\a\ 



\/3\ 2 = 1 : [Ai f ] = 11 



V 



Y(a,(3) 



1 (A^A) 

(At A) (At 2 A 2 ) 

F(a,/3)tM( 2 )(pg b) )F(a,/3), 
d \ 


^o(«,/3) 





/ 



^o(a,(3) 




aa 

a(3* 
13a* 
(3(3* 



(4.12) 



The dependence of A on a, (3 is left implicit. The superscript (2, 1) at the start of the above 
equations indicates that we are dealing with a general single-mode Mandel matrix obtained 
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from the two- mode Mandel matrix for the a — b system in the state p^\ by focussing on a 
variable linear combination A of a and b. As we will immediately see, for two-mode states 
both the 5 x 5 matrix M^(p^ b) ) and the 2 x 2 matrix M (2 '^(p^ b) ;a, (5) are important. 

The two-mode definitions of Mandel-type nonclassicality, sub-Poissonian statistics (sub- 
PS), super- Poissonian statistics (super-PS), etc are now as follows: 

{pg b) isQO-cl P(I a ,I b ,6)>0} 

=S> {M^\p ( f) > pifhas super-PS}; 
M (2) (pg 6) )^0 <=> {pg 6) is QO-noncl, has sub-PS}. (4.13) 

In the definition of super-PS, we used Eq. (4.11). The sub-PS case can be usefully separated 
into two types, depending on whether or not the nonpositivity of the 5x5 matrix M^(p^) 
is visible already at the single-mode level for some choice of coefficients a, f3. Thus we define : 

pg b) has Type I sub-PS M^ l \p^ b) ; a,/3)^0 for some a, # 

pg 6) has Type II sub-PS <=> M (2) (pg 6) ) £ 0, but 

MM^f'ja, p) > for all a, (3. (4.14) 

The physical meaning is that in Type I sub-PS, the Mandel level of QO nonclassicality is 
easy to detect already in terms of a suitable single-mode combination; while in Type II 
sub-PS, such nonclassicality is hidden or intrinsically two-mode in character [12]. 

For calculational purposes one can pass from the 5x5 Mandel matrix M^(p^) to a 
slightly simpler 4x4 matrix as follows. From Eq. (4.11), 



M (2) = Tr^H) 



J D 



1 Ct\ 
C B 

C = Tr(pg b) C0, B = Tr(p% b) B). (4.15) 




/ 



(When necessary the state will be indicated as argument of C,B). Then it is easy to see 
that 

M (2) (# } ) > <£> T = B-CC^>0, 

M (2) (# } ) ^ <£> r t 0. (4.16) 
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Thus the 4x4 hermitian matrix T determines whether we have super-PS or sub-PS. For 

the separation of the latter into Type I and Type II, we have for any complex 2-vector 

( jl' 



n 



0tM( 2 ' 1 )(# ) ;a,/3)0 = 

|0i + 2 CVo(«,/3)| 2 + |0 2 r^o(a,/3) t r^o(«,/3). (4.17) 
So given we are able to say : 

Type I sub-PS <£> /3) f r^o(a, P) < for some a, p; 

Type II sub-PS <£> fo(a, p^T^a, p) > for all a, p. (4.18) 

Indeed we easily find from Eqs. (4.12, 4.16) that 

det M^\p^; a, P) = if> (<*, P)^M<*, P)- (4-19) 

D. Examples of two-mode states and Mandel matrices 

We consider two examples. We have seen in Section II that a single-mode QO-noncl 
state, when combined with a second ancilla mode in vacuum (or in a coherent state) and 
then passed through a nontrivial U (2) BS, always results at the output in a two-mode state 
exhibiting NPT entanglement. We study this as the first example. 

The two- mode state in question is given in Eq. (2.16). It is understandable that its 5 x 5 
Mandel matrix is obtainable from the 2x2 Mandel matrix associated with the single mode 
input state p$ . Straightforward calculation shows that : 

= U(u){p^®\0) bb (0\}U(uy\ ueU(2): 

M( 2 )(pg b) ) = W( U yM^(p^)W(u), 

1 

U^Uu U^Un «21 M 21 

W{u)W( U y = l 2x2 . (4.20) 



W(u) 



Next using (4.12) we can immediately obtain the variable single-mode projection of this 
two-mode Mandel matrix : 

M^ 1 ) = Y(a,P^W{u)^M^\pi ] )W{u)Y(a,P) 

20 



1 



x lei 2 

f = Ull a + u 2 2f3- (4.21) 

From these expressions and the results of Section II, we find that the two-mode states pro- 
duced from single-mode states in the above manner have the following significant properties : 

(i) ff$ has QO-noncl PND =4> pffi has NPT entanglement; 

(ii) pj} has super-PS =4> fffi has super-PS; 
(Hi) p$ has sub-PS =>• pffi has Type I sub-PS, 

M^ip^-a^) t for every a, /3. (4.22) 

Of course only properties (ii) and (Hi) involve the Mandel matrix analysis; it is significant 
that in (Hi), every single-mode combination of the modes a and b displays sub-PS. To this we 
can add the following : states pffi obtained from states p$ via Eq. (4.20) can never display 
Type-II sub-PS; and any sub-PS in p$ leads to both Type-I sub-PS and NPT entanglement 

. ^ab) 

m p y D '. 

The second example is the two-mode generalization of the squeezed vacuum state defined 
for a single mode in Eq. (4.5). We take independent complex £, £' or u, u' and define: 

\^ ab \u,u')) = \^ {a) (u)) ® \^ {b) (oo')), (4.23) 

with the second factor involving an exponential in b^ 2 and b 2 applied to \0)b- This pure state 
is clearly also QO-noncl, but it is manifestly a product state of Schmidt rank one. Unlike 
the single mode case in Eq. (4.7), however, now the QO nonclassicality shows up at the 
Mandel level. The 5x5 Mandel matrix for the state (4.23) is easily found using Eqs. (4.6) 
and their 6-mode analogues : 



M«(#) 




M^(\^ ab \co,co'))) = 



1 

C B 



B = 



S 2 S' 2 

S 2 (2S 2 + C 2 ) 



S 2 S' 2 





S 2 S' 2 
e-^SCS'C 






e^SCS'C 
S 2 S' 2 




S 2 S' 2 ^ 



S' 2 (2S' 2 + C' 2 ) J 
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(4.24) 



Here S' and C are defined as in Eq. (4.6) but in terms of The 4x4 matrix T of Eq. 
(4.16) is: 



r = 



1 s\s 2 + c 2 ) 




V 





S 2 S' 2 e iT >SCS'C 
e-^SCS'C S 2 S' 2 

S' 2 (S' 2 + C' 2 ) 







->/ 2 / of 2 



(4.25) 



The eigenvalues of T are S 2 (S 2 + C 2 ), S' 2 (S' 2 + C' 2 ), SS'(SS' + CC) and SS'^S" - CC"). 
Assuming that £, £' are both non zero, the last eigenvalue is negative, leading by Eq. (4.16) 
to the conclusion that M (2 \\^ (ab \uj, u'))) ~£ or that the state \^ ab \uj, u')) has sub-PS. 
This is an interesting and somewhat nonintuitive result since we have seen in Eq. (4.7) that 
each factor in the product state \^ ab \u,uj')} has super-PS. We must now see whether it is 
Type I or Type II. For this we must compute the 'expectation value' of T in Eq. (4.25) for 
the four-component column vector ipoi®, P) as required by Eq. (4.18) : 

^o(tt,/3) t r^o(tt,/3) = 

(\a\ 2 S 2 + \/3\ 2 S' 2 ) 2 + \a\ 4 S 2 C 2 + \/3\ 4 S' 2 C' 2 + 2SCS'C'^{e ir >{a* (5) 2 ) 

> (\a\ 2 S 2 + \(3\ 2 S' 2 ) 2 + (\a\ 2 SC - \[5\ 2 S'C') 2 > 0, (4.26) 

since $l(e ir > (a* (3) 2 ) > -\a\ 2 \/3\ 2 . It follows that the sub-PS of the product state \^ ab \co, u')} 
is of Type II, it is hidden or intrinsic. This is consistent with the fact that the individual 
factors and \^ b \u')) are both super-PS. 



V. TWO-MODE MANDEL LEVEL NONCLASSICALITY TO ENTANGLEMENT 
BY BS ACTION 

We now consider a two-mode state fffi which is QO nonclassical and of such a nature 
that this property is seen at the level of the Mandel matrix, i.e., M^(p^) ^ 0. In such a 
case, even if is of product or separable form, its passage through a U (2) BS could result 
in an entangled state, possibly of NPT type. Our aim is now to see quantitatively how much 
beyond the nonpositivity of the Mandel matrix one has to go to reach NPT entanglement. 
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We set up the general framework to examine this, then illustrate it by an example. For 
simplicity we use a 50:50 BS rather than one corresponding to a general u e U(2). 

We choose the U(2) element and corresponding unitary operator action as follows [11, 34] : 



1 1 



u 



e U{2) : U, 1 



^ a a) \ ina 
b W \b 6t 



V2 \ _i i 

At the operator level, action by conjugation on C, C\ B of Eqs. (4.8, 4.11) is: 



U^CUo = V C, U^&Uq = C% T , U^BUq = V BV T , 



Then when the state pffi is transformed by this BS action to 



( 1 1 11^ 
-1 1 -11 
-1-1 1 1 
1 -1-11 



V 



-/ (a6) _ f j 4ab)fj-l 
P D — U 0Pd U > 



the new Mandel matrix is given by a transformation using V$ : 

1 C7+ 



M^pjf) 

M (2) (/3 ,(a6) ) 




/ (aft) 



1 # 



'i o N 

v° V ° T ; 



C = V C, B' = VqBVq . 



Thus T' is related to Y by congruence : 



Y' = B'- C'C* = V YV*. 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



To test next whether p' ^ is NPT entangled, we pass to its partial transpose p' ^ PT 
and evaluate the 'expectation value' of a suitably chosen nonnegative hermitian operator 
with respect to it. If this turns out to be negative, then the output state p'ffl is definitely 
NPT entangled. To construct such a test which involves as closely as possible the use of 
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M ( 2 )(p'g 6) ), hence of M^ 2 \p^ b) ), plus something additional, we should use a 'matrix of 

( i Wi ct) 

operators' similar in structure to v 7 , i.e, making up a 'hermitian nonnegative' 

matrix of operator entries, such that when the partial transpose operation is switched from 
p (ab)PT ^ o ^ s 'matrix', we obtain the expectation values of C, C*t and B in p >( £ b \ plus 
something additional. Now we have seen in the passage from Eq. (3.7) to Eq. (3.8) that 
the PT operation converts b^b k to tf k b> , and b>w k to b k b^ . Keeping these motivations and 
facts in mind we construct a 5 x 5 matrix of operators as follows : 



E 




db 
1 St 



it 



PT 



Y 



ata db aW b'b 



1 ct\ 

C £ J 


+ 

V 


\ 
Y ) 




a f a 




















S+a 


ata 


+ fetfe H 


- 1 ot a 







atfe 


&t& 



(5.6) 



We see that in the process of expressing the various operators involved in normal ordered 
form, as anticipated an additional piece Y linear in the entries of C appears. Then a test 
for NPT entanglement of p' ffi is to evaluate 



rp / ~/ (ab)PT 

Tr(p K D 



1 Bt 



Tr(p'Sf 



E 



E 



i it 



PT 
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\ 










Y' = 



C[ 




c> 2 










V 



c 3 





Ci + Ci + l c 3 




c 2 



C 4 j 



(5.7) 



and see if this matrix is indefinite. By Eq. (5.4), the complete 5x5 matrix here is a 
congruence transformation applied to the initial state Mandel matrix M^(p^) plus a 
4x4 piece coming from Y, namely it is : 




V C V BV T + Y' 



(5.8) 



Therefore by Eq. (4.16) the positivity or otherwise of the matrix (5.7) is equivalent to the 
positivity or otherwise of either of the two following 4x4 matrices at the level of T : 



rT 



T + V T Y'V , 



V nVJ = VoTV^ + Y'. 



(5.9) 



Nonpositivity of either Q or VqQVq is proof of NPT entanglement of the output state 
p £) b \ It is interesting to see the precise quantitative manner in which this test goes beyond 
examination of M^ 2 \p^) or V alone. 



A. An illustrative example 

To see how the general procedure developed above works, we study a family of states 
which is analytically quite simple. We begin with the family of two-mode pure states of 
infinite Schmidt rank, 

co n 

lfl) = e -^\ 2 J2-^=\n,n), peC, (5.10) 
n=o V"! 
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form the density matrix p( ab "> = \p}(p\, and pass to pffi via Eq. (3.3) : 



= e~ x J2 —r \n,n)(n,n\, A = |/i| 2 > 0. 



?1=0 



(5.11) 



This is clearly separable though not of product form. For the Mandel matrix analysis, | /x) {/x | 
and fffi are equivalent. 

The matrices C, C\ B, F involved in M^(p^) are easy to calculate since 



(ata> = (S+S) = A, (a t2 a 2 ) = (bH 2 ) = A 2 . 



(5.12) 



We thus have : 



C 








/ A 2 



B 



A 2 + A 



A 2 + A 



A 2 + A 



A 2 + A 



A 2 



OA 

A 2 + A 

A 2 + A 

AO 



(5.13) 



The eigenvalues of T being A(A + 1), A(A + 1), A, —A, the state in (5.11) is QO-noncl. 
To find its type we compute 



Vo(«, P)^ (a, (3) = 2|a| 2 |/3| 2 A(A + 2) > 0, 



(5.14) 



so these states display hidden or Type II sub-PS. 

In passing we note that the state p^ of mode a obtained from Eq. (5.11) by tracing over 
b alone is 



4ab) _ —A 



A T 



n=0 



(5.15) 



for which the diagonal weight P(I a ) is 



P(I a ) = 6(I a -X). 



(5.16) 



Partial trace over a gives exactly similar results for mode b. Thus both p$ and p$ are 
QO-cl, with their PND's coinciding exactly with that of a coherent state. 
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Now we pass the two-mode state pf^ of Eq. (5.11) through the BS U of Eq. (5.1); the 
resulting p' ^ is 



P D 



= e 



£ 

n=0 



a 



i t2 -6 t2 ) n |0,0)(0,0|(a 2 -6 2 ) n . 



(5.17) 



To apply the NPT entanglement test based on Eq. (5.9) it is convenient to examine VqVLVq . 
Combining Eqs. (5.4, 5.13) we find the matrices T', Y' associated with p'ffl to be 



r' = v rv^ 



'±A 2 +A 







-§A 2 \ 









-\x 2 -x 




A 2 - 


-A 


±A 2 


{ -\* 







o \x 2 + x J 


f X 


°1 
















2A + 1 









v o 


A J 







Y' 



Therefore according to Eq. (5.9) we have to test the positivity or otherwise of 



(5.18) 



V QV^ = VqTV^ + Y' 
( ±A 2 + 2A 



V 



o 
o 





|A 2 
-±A 2 - 






-|A 2 -A 
A ±A 2 + 2A + 1 




±A 2 






|A 2 + 2A 

1 \2 



(5.19) 



The (2,3) submatrix here is indefinite as it has determinant — ^A . This establishes that 
p'n^ °f Eq. (5.17) is NPT entangled. Qualitatively speaking, even though Y' in Eq. (5.18) 
is nonnegative, the total matrix VqQVq in Eq. (5.19) is indefinite, with V dominating Y' . 
The emphasis here has been to show that the NPT entanglement produced by BS action can 
indeed be witnessed by the test based on Eq. (5.9), which goes beyond the examination of the 
Mandel matrix in a precise manner. 

We can characterize the NPT entanglement we have proved in this example further. From 
the expression in Eq. (5.17), the terms for n = and n = 1 are respectively: 



' 10,0)(0,0|, 



- A (|2,0>-|0,2»((2,0|-(0,2|), 
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(5.20) 



giving the matrix elements 



\P d Joo,oo — e , 

[P D J20,20 \P D J02,02 — ~ \P D J20,02 — — IP D J02,20 



A 



—A 



(5.21) 



One also obtains from the n = 2 term in Eq. (5.17) the matrix element 

A 2 



(n' ( a0> \ 
[P D J22,22 



—A 



(5.22) 



To demonstrate distillability we follow the recipe of Ref. [40], and project p'ffl into the 
2x2 bipartite subspace spanned by |0), |2) of the a- mode and |0), |2) for the 6-mode. The 
resulting 2x2 state written in the basis |0,0), |0,2), |2, 0), |2, 2) reads 



(ab) 
P D ^ 



2x2 



^ 1 ^ 

o I -I o 

o -I I o 

o o o 4^ 



(5.23) 



That (yf)'^ ) 2 2 is inseparable is manifest, for its partial transpose has negative determi- 
nant. Indeed, the partial transpose has an eigenvector of the form a|0,0) + /3|2,2) with a 
negative eigenvalue. This demonstrates that the NPT entanglement of p'ffl in Eq. (5.17) 
is in fact distillable. 



B. The two-mode squeezed vacuum 

This state defined in Eq. (4.23) is both pure and of product form, i.e., of Schmidt rank 
one. It was studied in the previous Section via its Mandel matrix, which showed it to be 
QO-noncl with Type-II (i.e., hidden) sub-PS. In principle, after passing this state through 
the BS U , we can apply the test based on Eq. (5.9) to see if the output state is NPT 
entangled. However, this involves some amount of algebra. Fortunately, the fact that the 
output state of the BS is entangled can be seen in this case simply by inspection and without 
any calculations: 

W (a6) (^')> = t>oexp{i(£at 2 - Ca 2 ) + \(^ 2 - tb 2 )}^^) 

= exp{i(£(at - ftt) 2 - C(a - bf) + i(£'(at + Pf - f'*(a + 6) 2 )}|0, 0). (5.24) 
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This is because the final unitary operator acting on 1 , 0) is clearly not the tensor product 
of individual unitary operators acting separately on the two vacua. It is of course important 
that at least one of the factors \ip^(u)), \^ b \u')} in the initial product (4.23) be QO- 
noncl. A two-mode pure product QO-cl state is necessarily a product of single-mode coherent 
states, and this product structure is maintained by BS. 



VI. TWO-MODE NONCLASSICALITY TO THREE-MODE ENTANGLEMENT 

In the preceding Section we studied the possibility of a U (2) BS converting a two-mode 
QO-noncl separable state into an entangled one since both nonclassicality and entanglement 
are meaningful concepts for such systems. The entire discussion was within the framework 
of the space of states of a two-mode system. 

Now we present a treatment of two-mode states analogous to that in Section II for single- 
mode systems. That is, we couple a given two-mode state pffi to a third ancilla mode 
in vacuum, pass such an input state through a l U(3) beamsplitter' (a classicality 

preserving passive system), and obtain a three- mode output state Pout^- We then test 
whether this shows NPT entanglement as a consequence of (Mandel level) QO nonclassicality 
assumed to be present initially in p^; the PT operation is applied to the c-mode. The 
motivation is to explore the algebraic expressions and form of the test one is led to, apart 
from carrying the physical process described in Section II to the next higher level. 

We begin with pffi for which M^(p^) shows QO nonclassicality. With the ancilla 
c-mode in vacuum we have an input three-mode state 

pff c) =P&%|0> CC (0|, (6.1) 

strictly analogous to Eq. (2.15). To a general matrix u G £7(3) we associate a passive 
'beamsplitter' which unitarily mixes the annihilation operators of the three modes in a 
manner analogous to Eq. (2.13), now conserving N a + N b + N c . In the three- mode Hilbert 
space this BS u acts through a unitary operator U, and we have [19, 20] 



u = 



( u n un u 13 \ 

«21 «22 «23 
\ «31 «32 U33 ) 



e 17(3) -> U : U ] U = UU^ = 1, 
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u 



b 



U 



-1 



u 



U 



u- 



u = 



u 



U(N a + N b + N c ) 



b 



, u- 1 



/at\ 

St 
Vet; 

/at\ 

St 

yet; 



U- 1 = u 



U = u* 



/at\ 

St 
\& j 

/at\ 

St 
Vet; 



(N a + N b + N C )U. 

Therefore upon passage through this BS the state in Eq. (6.1) changes to 

pff = Up^U- 1 = UVp^ ® |0) CC (0|}L>- 1 . 



(6.2) 



(6.3) 



To test this output state for NPT entanglement, we apply the PT operation to the c- 
mode and then evaluate the 'expectation value' of a suitably chosen hermitian nonnegative 
operator : 

A = a + a\hc + aj)c + a^a}& + ajflct : 
Tr(pit c)PT A^A) = a^Xa, 
Tr(p^ PT ( 1 N 



X 



S+ct 
ac 



1 ac Sc at c t Stgt ^); 

cci 

m = « 2 . (6.4) 
a 3 
\ « 4 ; 

In spirit the approach here is similar to that used in the preceding Section leading to Eq. 
(5.7), though here we are dealing with the extension from two to three modes. The 5 x 
5 hermitian matrix X, constructed by taking entry wise expectation values as defined, is 
designed to be related to the input Mandel matrix M^(p^) but going beyond it in a well 
defined way. Developing it we find : 

PT 

1 a< be -rc : S T /-t ^ 

> = 





( 1 \ 




atct 


< 


S+ct 




ac 




V hd J 
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1 1 N 

tic 



1 &a &b a)c tic 



+ 



V 





















z = 



a) a + $c + 1 tia 

a)b ttb + &c + 1 

+ (1 + c f c 



1 
1 



(6.5) 



Use this in Eq. (6.4), implement the conjugation U 1 {- ■ -)U and use the fact that the c-mode 
is initially in vacuum to get : 

/ 1 \ 



X 



(ab) . 



a'tc' 



( 



1 c'ta' a'tc' fo'tc' ) : 



(ab) 



















Z' 



), 



Z' 




t St 



+ (i + etc') 



1 
1 



(a'\ 

V 



U 2 1 U 2 2 
\ u 31 U 32 J 



(6.6) 



The appearance of the extra terms Z, Z' is a result of normal ordering similar to the 
appearance of Y in Eq. (5.6). One can now disentangle the u- dependences and express the 
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result in terms of M^(p^) and an additional piece involving C = Trip^C) : 



X = W{u)M( 2 \p<g b) )W{ U y + 



W(u) 












: 



o o 




Z' 








u 



11 



u 31 U 32 



u 



12 



u 31 ^32 



21 



U 



22 



Mil Wl2 



U21 U 2 2 



z' 
z' 




u 



32 



Ull U12 



U21 U 2 2 



at fet 



^ u *l M 21 ^ 



V M 12 M 22 / 
/ 

+ (1 + (u^ 1 a t + w^Xwgia + m 32 6)) 



1 
1 



(6.7) 



If this matrix X, dependent on pffi and u £ U (3), is indefinite, Pout^ °f Eq. (6.3) is NPT 
entangled. For this to happen, as we have assumed, pffi must be QO-noncl, since the BS 
U would map any QO-cl input into similar output, the ancilla being in a QO-classical state 
(vacuum) . 

A. Two illustrative examples 

The first is a two-mode state with only a finite number of photons, so that its QO 
nonclassicality is a foregone conclusion: 

Jab) 



p% 0> = p|2,0)(2,0| + g|l,l)(l,l|+r|0,2)(0,2|, 
P, q,r>0 p + q + r — 1. 



(6.8) 
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This is separable, though not a product state. The only non vanishing expectation values 
needed to construct the Mandel matrix are 



(a f a) = 2p + q, (Pb) = 2r + q, (a f2 a 2 ) = 
(tfattb) = q, {P 2 b 2 ) = 2r. 

Therefore the Mandel matrix is 

1 q + 2p q + 2r 

q + 2p 2p q 

q 

q 



2p, 



(6.9) 



(6.10) 



\q + 2r q 2r 
The determinants of various nontrivial 2x2 submatrices, the one nontrivial 3x3 submatrix, 
and finally of Af( 2 \pffl) itself, are (indicating the submatrices by the relevant rows and 
columns) : 

(1,2): 2p-(q + 2p) 2 ; (1,5): 2r - (q + 2r) 2 ; (2,5): Apr - q 2 ; 

(1,2,5): q 2 -4pr; detM (2) (p^ 6) ) = q 2 (q 2 -4pr). (6.11) 

One can easily visualize situations for which the (1,2) and (1,5) submatrices become in- 
definite, for instance q close to unity and p,r close to zero. In any case, since the (2,5) 
sub determinant is opposite to the (1, 2, 5) subdeterminant in sign and also to the full deter- 
minant, the state in Eq. (6.8) is always QO-noncl at the Mandel matrix level. 

The type of sub-PS can be easily determined. From Eq. (6.10) we find the 4x4 matrix 
T to be : 



s„ 



Therefore also 



q- (q + 2p)(q + 2r 

q 

q 

\q-(q + 2p)(q + 2r) 5 b 

5 a = 2p - (q + 2p) 2 , 5 b = 2r - (q + 2r) 2 . 



M<*,PyrM<x,P) = 2p\a\ 4 + 4q\a\ 2 \f3\ 2 + 2r\ 



\ 



(6.12) 



-((q + 2p)\a\ 2 + (q + 2r)\ 



|2\2 



(6.13) 
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For a — 1, /3 — this is S a ; for a = 0, /3 = 1 it is 5&. We now consider p running over its 
range [0, 1] in successive portions and draw corresponding conclusions: 

p = : q = 5 b = -2; q>0 =>■ 5 a < 0; 
$<p<\: S a >0 2p- (p-r + 1) 2 > 

(p - rf + 1 - 2r < =4> 2r > 1 5 6 < 0; 
S a = (p - r) 2 + 1 - 2r = p ^ r, 

2r > 1 =>■ 5 6 < 0; 

1 1 
p=-: g = =^ p = r = -, 5 a = 5 b = 0; 

g > 5 a <0; 

- <p < 1 : 2p > 1 5 a <0. (6.14) 
2 

Thus in every situation except p — r — 1, q — 0, either 5 a or 5 b is negative. In this one 
exceptional case we find from Eq. (6.13) : 

p = r = ^q = 0: ^ (a, ^^(a, /3) = -2|a| 2 |/3| 2 , (6.15) 

which is negative for a, (5 ^ 0. This establishes that the state (6.8) is of Type I sub-PS. 

Now we couple this state to the third c-mode in vacuum, and pass it through a particular 
U(3) BS, namely a 50 : 50 BS acting on the b and c modes alone. The output state is 
calculated using Eq. (6.3), and to test whether it is NPT entangled we need to calculate 
the matrix X of Eq. (6.7) involving the Mandel matrix term and the added Z' term. The 
choice of u e U(3), the resulting W(u), and the two parts of X are as follows : 





(l 


^ 




u = 





l/y/2 l/y/2 


e 17(3); 




1° 


-l/y/2 l/y/2) 






(l 





\ 







-l/y/2 





W{u) = 








-1/2 







-l/y/2 




1° 





-1/2) 
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W(u)M {2 \p^ b) )W(u)^ = 



1 





z' = 



V 



o 






-r - q/2 





-r - q/2 


q/2 








u 


n 


r/z 


n 

u 


r/9 


u 





q/2 


n 
u 





r/2 





r/2 


_l_ r _ 


f 1 








q + 2r + 1 J 





(6.16) 



The dotted lines indicate where the 2x2 block Z' has to be inserted. Leaving out the trivial 
second and fourth rows and columns as they do not couple to any others, the determinants 
of the various 2x2 submatrices and the 3x3 submatrix in X are : 



(1,3): r/2- (g/2 + r) 2 ; (1,5): 5r/2 + q + 1 - (q/2 + r) 2 ; 
(3,5): r(g + 2r + l)/2; (1,3,5): (q + 2r + l)(r/2 - (q/2 + rf) 



(6.17) 



Comparing these with Eqs. (6.11) we see: whenever the QO nonclassicality of pffl manifests 
itself in the (1,5) submatrix of M^(p^) being indefinite, simultaneously the 3-mode state 
Pout^ displays NPT entanglement. If on the other hand the (1,2) submatrix of M^(p^) 
were indefinite, then by suitably altering the U(3) element u in Eq. (6.16) we can again 
achieve NPT entanglement of Pout^- 111 either event, we see how a U (3) BS can produce NPT 
entanglement starting from a two-mode nonclassical state (6.8), and how the signatures go 
beyond the indefiniteness of M^ 2S> (p^) in a precise manner. 

The second example to illustrate the ideas of this Section is similar in structure to the 
example (5.11) of the preceding Section, but differs in certain details. For a real nonnegative 
parameter 77 we define the separable state 

J2n 



Pd 



1 00 

-Y 
c z - 



\n, n)(n, n\, 



(6.18) 



n=o (2n)! 

where C = Cosh i], and later S = Sinhr/ and t = tanhr/. The case i] = corresponds to the 
two-mode vacuum, and so we take < i] < 00. Using the elementary sums 

n 2n 



]T(n or n 2 )-^^ = 'jS or 'j(S + r]C), 



rj 



(6.19) 



n=0 
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the nonzero expectation values needed for the Mandel matrix are : 



(at«) = (StS) = |f; 



(6.20) 



Therefore we find : 



rjt 
2 



f Uv-t) 



















2 



\ 



(6.21) 







o *(v + t) 
2 4 V* ■ v o o 

Leaving out the third and fourth rows and columns, the remaining 2x2 subdeterminants 
are : 



Uv + t) 






* 2 v7 + t) 



;i,2)and(l,5): ^ - 1); (2,5): 



(6.22) 



The function ^ — t decreases monotonically from to —1 as r] runs from zero to infinity. 
We see that the state (6.18) is QO-noncl for all rj > 0. To determine its Type we compute 
T and its 'expectation value' in ip (a,(3) : 



C 2 



4 



-t 







-9- 4- 1 

c 2 ^ 1 



\ 



rj + t 
rj + t 

n ji / 



y ^ + t 



M<x,P)*rM<x,P) = |{^-t + 2|a| 2 |/3| 2 (r/ + 3t)}. 



(6.23) 



At both a = 1, /3 = and ct = 0, (5 = 1 the last expression is negative, so the state (6.18) is 
QO-noncl Type I sub-PS. In this context we note that the single-mode state obtained 
from (6.18) by tracing over b is 



C^(2n)! 

and this has the Mandel matrix and determinant 



\n) aa (n\, 



(6.24) 



M (1) (#) 



1 nl 

2 

rjt rj 



detM«(#) = If(JL- t )<0. 



77 77 



4 V C 2 



(6.25) 
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The properties of p$ are identical. So in contrast to the state (5.11), now both p$ and p^ 
are QO-noncl, accompanying the Type I nature of pffi . 

We now apply the NPT entanglement test outlined in Eqs. (6.4, 6.6, 6.7). The necessary 
expressions are : 



W{u)M^(p^)W( U y 



1 





vt 
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nt 
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Mv- 
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Mv- 


















v t 
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Mv- 
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Z' 



1 + 



3ijt 
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o 



l + 25 / 



(6.26) 



The 2x2 matrix Z', which is positive definite, has to be 'added' at the lower right hand 
corner of the 5x5 matrix, leading to X of Eq. (6.7). Then the positivity or otherwise of 
X has to be examined. However, even without taking account of Z', and unaffected by Z', 
the (1,3) subdeterminant of X is ^(qi — t), which is negative. This establishes the NPT 
entanglement of Pout^ i n this example. 

The considerations of this Section show that the scheme described in Section II, elevating 
single mode QO-noncl states to the two-mode level and then allowing BS action to produce 
NPT entanglement, generalizes to the next higher level. The method of Mandel matrices is 
a practical way to see these processes in action. 



VII. GENUINE TRIPARTITE ENTANGLEMENT 

Now that the main methods of our approach-signatures of nonclassicality at the Mandel 
matrix level, their nontrivial extensions to signatures of (NPT) entanglement created by 
BS action - have been applied to several examples, we go on to consider some more subtle 
features of entanglement. We will show via an example that in the three- mode case the 
BS action on an initial two-mode nonclassical state can lead to genuine residual tripartite 
entanglement. This is in the sense of [43], whereby the output is a tripartite state similar 
to the GHZ state [44] : there is no bipartite entanglement when any one of the three modes 
is traced away. 
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We consider the state (5.11) studied in Section V, and subject it to the treatment of 
Section VI. As we have seen, this (separable) state shows Type II sub-PS. With pffi as in 
(5.11), we pass the state (6.1), 



p ( D b) ® |0)cc(0| = e~ A ]T ^-\ n ,n} abab (n,n\ <g> |0) cc (0|, 



A 



n=0 



n 



(7.1) 



through a 50:50 b — c BS, a special Z7 (3) element, whose action on the mode operators b and 
c is 



(7.2) 




The resulting state is 



pZ c) = uip™ ® |o) cc (o|)f/- 1 

oo 



(n\®(P + $) n \0,0) bcbc (0,0\(b + c) n 



ra=0 

00 A n n! 



= e 



■E 

n=0 



yn 



\n) nn \n\ 



^ |r,n - r) bcbc (s,n- s\ 



(7.3) 



r ,s=o yr\(n — r)!s!(n — s)! 

Clearly this is separable in the a/bc cut. However it is entangled in both the c/ab and b/ac 
cuts as we show below. As a test for NPT entanglement in the c/ab cut, we evaluate the 
'expectation value' of a suitably chosen positive operator on the partially transposed output 
P^ont >PT ■, the partial transpose being applied on the c mode. For the choice of operator A^A 
where 



A = «o + ctibc + a^a^a, 
a test for NPT entanglement would be to check for violation of positivity of 
Tr(pff PT iti) = Tr(/ff (iU)^) 

( a* a\ a* 2 ) X 



(7.4) 



u \ 



a 1 
\a 2 j 
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X = Tr(p, 



(abc) 
out 



feet 



fotg b^bc^c b^ca^a 
at a a^abc^ a^aa^a 



(7.5) 
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Using the fact that initially the c-mode is in the vacuum, we find : 

' 1 A/2 A \ 

X= A/2 A 2 /4 A(A + l)/2 . (7.6) 
v A A(A + l)/2 A(A + 1) j 

As the (2,3) submatrix of X has negative determinant, Pout^ i s NPT entangled across the 
c/ab cut. It is easy to see that a similar test with the same choice of A, except that now the 
PT operation is applied to the b mode, yields the conclusion that piut^ i s NPT entangled 
across the b/ac cut. So we find in this example bipartite entanglement in a tripartite setup, 
as a result of BS action. 

Now to show that the entanglement is genuine tripartite, 'residual' in the sense of [43], 
we have the following : 



P-t = Tr c (p^ ut >) = e ^ \ n ) aa (n\ ^ 

n=0 Z r=0''V i ' )■ 

Pout = Tr 6 (p^ ut ; ) = e ^ -^r\n)a a(n\ 2^ 



n=0 



2" 1 r tS r!(«- r)P 



pS = ^off) = «-» £ ^ £ |r ;"~ r) ^ < ' , ""» - (77) 

n =o z r,s=o yr!(n — r)!s!(n — s)! 

The first two are manifestly separable. It may not be obvious at first glance that the third 
is also separable but a closer look shows that it can be written in the form 

PoS = e- A t>(£ ^\n) bb (n\ ® 1 > cc <0 1 ) 1 , (7.8) 

where C/ is the 50:50 b — c BS (7.2). Note that J2™=o ^\ n )bb( n \, the state of the b- mode at 
the input, is simply the phase averaged (coarse grained) version of the coherent state \ V\}- 
Thus p^t is the outcome of a classical state passed through a BS, so it is classical and hence 
separable. 

It is interesting that the feature of genuine tripartite entanglement is reminiscent of Type 
II sub-PS for a two-mode state at the Mandel level, where the nonclassicality never shows 
up at any single mode level. An interesting question in this context is the possibility of 
extension of monogamy relations to this non- Gaussian case [45, 46]. 
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VIII. FURTHER PROPERTIES OF MANDEL PARAMETERS AND BEAM- 
SPLITTERS 



For one and two mode field states, we have used the 2x2 and 5x5 Mandel matrices 
respectively to classify the states in a physically useful manner. It is convenient to also have 
suitably normalized single parameter - 'scalar' - measures of nonclassicality defined in terms 
of the Mandel matrices. In the two-mode case a useful requirement would be invariance of 
such measures under BS action. 

We begin with the single mode case. Here the Mandel Q parameter was defined in [8] as 

(Ai\y 2 - (N a ) 



Q 



(N a ) 
(qt 2 q 2 ) - (qtq) 2 

detAfWGSj?) 



(8.1) 



Here all expectation values are in the state p$, and Eq. (4.2) has been used. From the 
classification (4.4), Q > and Q < correspond respectively to super-PS and nonclassical 
sub-PS cases. This parameter is bounded below by —1, which is a convenient normalization. 
For Q > there is no upper bound. 

In attempting to generalize to two modes, as a first step we show that the separation of 
QO-noncl states into Types I and II is BS action invariant. For any u G U(2), from Eq. 
(2.13) and the direct product structure of C in Eq. (4.8) we easily obtain: 

(j-^CU = VC, u-^u = &V\ 
U~ X BU = VBV\ 

V = u*®u. (8.2) 

Therefore, with C and B defined as in Eq. (4.15), under general BS action we have: 

Pd 



~> 6 ) = UpWjj- 1 ^C' = VC, B' = VBV\ 



r' = VTV ] . (8.3) 

Now in the QO-noncl family of states, corresponding to M^(p^) ^ or equivalently 
to T ^ 0, the further separation into Types I and II is given in Eq. (4.18). Here it is 
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the 'expectation values' of T in four component column vectors ipo(a, (3) that are relevant. 
However these vectors too have a direct product structure (4.12), so they are mapped into 
similar vectors under the above changes : 












T 

= u 











(8.4) 



This proves that the separation into Types I and II is preserved under BS action. More 
specifically, if a BS converts a QO-noncl separable state of a definite Type into an NPT 
entangled state, this change occurs within the subfamily of that Type. 

Now we generalize (8.1) to the two mode case. Keeping the requirement of BS action 
invariance in mind, we define the two-mode Mandel parameter as 

Tr(r) - ||r|| 



Q' 



.5) 



2«ata) + (StS))' 

Here ||r|| is the trace norm of T, which for hermitian T is the sum of the absolute values of 
its eigenvalues. Therefore Q' is simply the sum of the negative eigenvalues of T divided by 
the expectation value of the total number operator N a + iV&. From Eq. (8.3), the invariance 
of Q' under BS action is obvious. 

The parameter Q' vanishes for QO-cl states (as defined via the Mandel matrix), and 
is strictly negative for QO-noncl states. For two-mode product Fock states, for instance, 
Q' = —1. To follow the distinction between the two Types, we combine Eqs. (4.12, 8.1) to 
define the variable single mode Mandel parameter 

(A^ 2 A 2 ) - (A^A) 2 



Q(a,0) 

A 



aa + (5b. 



The minimum value of Q(a, f3) as a, [5 vary is also useful : 



Q min = min Q(a,P). 



(8.6) 



From all the previous discussions we draw up a table of results characterizing various two- 
mode states (always at the Mandel level) in Table I. For some of these, we have examples 
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Category 


Definition 


Description 


(i) 


Q min > 0, Q' = 


QO-cl 


(iia) 


Qmin < Q, < q 


QO-noncl Type I 


(iib) 


Q' < Q min < 


QO-noncl Type I 


(iii) 


g min > o, q' < o 


QO-noncl Type II 



TABLE I: 

from previous Sections. All QO-cl states, including the states (4.20) for detM( 1 )(pg ) ) > 0, 
come under category (i). On the other hand, the states (4.20) for detM^(p^) < belong 
to category (ii). The subclassification into (iia) and (iib) is subtle, but simple examples of 
each can be provided. For category (iia) we consider the product of a Fock state at the a 
mode and a coherent state at the b mode : 



= \n) a <8> \z) b . 

This QO-noncl separable. The matrix T is diagonal, 



(8.8) 



T = diag(— n, n\z\ 2 , n\z\ 2 , 0), 



J.9) 



leading to 



Q 2 = -n/(n+ \z\ 2 ) > -1. 



(8.10) 



On the other hand, Q(1,0) = —1, so <5 min < Q' < which falls under (iia). For category 
(iib) we can take the states (6.18) which are QO-noncl and separable. Using Eqs. (6.20, 
6.23) we find : 

Q' = -1/2 : 



= Y t {^-t + 2\a\ 2 \(3\ 2 ( V + 3t)}. 
The minimum of Q(a, (3) is reached when either a or /3 is zero : 

o min = —\^L n > -I 

As > Q min > Q', this falls under category (iib). 
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(8.11) 



(8.12) 



Another interesting example for (iib) is the class of pure states obtained as an equal 
in-phase superposition of product Fock states with given total occupation number n : 

1 n 

\^) = ^—T.\^-r). (8.13) 

For the cases n — 1,2,3,4 the numerically computed values of Q' are —1, —1.085, —1.123, 
— 1.143 respectively. Since in any case Q(a,f3) and Q mm are bounded below by —1, we have 
Q' < Q mm . We may also note that these entangled states cannot be produced by BS action 
on product Fock states, except when n — 1. 

Turning finally to category (iii), we have examples from Sections 4 and 5. For the two- 
mode squeezed vacuum state (4.23), we have using Eq. (4.26) and the eigenvalue spectrum 
of T stated after Eq. (4.25) : 

Q' = SS'{SS' - CC')/{S 2 + S' 2 ) < 0, 
Q min > 0. (8.14) 

For the family of states (5.11), from Eqs. (5.12, 5.13) we find: 

Q' = -l/2, Q min >0. (8.15) 

So in both cases we have category (iii) states. 

IX. CONCLUSIONS 

In this work we have investigated the relationships between quantum optical nonclassical- 
ity of the phase insensitive type and entanglement in multimode radiation fields In particular 
we have examined the possibilities of converting nonclassicality in such fields into entangle- 
ment through the use of classicality preserving passive devices such as beamsplitters. For 
the case of a single mode, after giving a complete characterisation of the quantum optical 
nonclassicality at the level of phase insensitive quantities, we have shown that such states 
through a beamsplitter action with vacuum or more generally a coherent state at the other 
port, always give rise to an NPT entangled state. For the case of two mode radiation fields we 
have presented a test which simultaneously witnesses both nonclassicality and entanglement 
in such states and have also developed a scheme based on Mandel matrices for characterising 
and classifying nonclassicality in one and two mode states. In particular, it is shown that 
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in the two mode case, the characterisation at the level of Mandel matrices permits us to 
divide two mode states into two categories-Type I where the nonclassicality manifests itself 
already at the single mode level and Type II where the nonclassicality is intrinsically two 
mode in character with no Mandel type signatures of nonclassicality at the single mode level. 
We have also examined in detail the possibility of a Z7" (2) beamsplitter converting Mandel 
level nonclassicality in a two mode separable state into an NPT entangled state and have 
given tests for NPT entanglement in the state resulting from such an action. Distillability 
of the state so produced is demonstrated in one case. Further, in a similar spirit as for the 
case of a single mode, we have also analysed the action of a U(3) beamsplitter on a two 
mode Mandel level nonclassical state and have derived conditions under which the two mode 
nonclassicality manifests itself in NPT entanglement in the resulting three mode states. In 
this context, we have also shown, via an example, how such an action can lead to a genuine 
tripartite entangled state in the sense that there is no bipartite entanglement when any one 
of the three modes is traced away. With a view to ease in categorisation of nonclassicality 
in two mode states, by appealing to invariance under beamsplitter action we have suggested 
analogues of the Mandel Q-parameter originally introduced in the context of single mode 
radiation states. 
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